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We analyze the effects induced by the bulk (or second) viscosity on the dynamics associated to 
the extreme gravitational collapse. Aim of the work is to investigate whether the presence of viscous 
corrections to the evolution of a collapsing gas cloud influence the top-down fragmentation process. 
To this end, we generalize the approach presented in the Hunter work [l[ to include in the dynamics 
of the (uniform and spherically symmetric) cloud the negative pressure contribution associated to 
the bulk viscosity phenomenology. Within the framework of a Newtonian approach (whose range 
of validity is outlined), we extend to the viscous case either the Lagrangian, either the Eulerian 
motion of the system addressed in [I[ and we treat the asymptotic evolution in correspondence to a 
viscosity coefRcient of the form C = Co P^^^ [p being the cloud density and = const.). 

We show how the adiabatic-like behavior of the gas is deeply influenced by viscous correction when 
its collapse reaches the extreme regime toward the singularity. In fact the density contrast associated 
to a given scale of the fragmentation process acquires, asymptotically, a vanishing behavior which 
prevents the formation of sub-structures. Since in the non- viscous case the density contrasts remain 
constant, we can conclude that in the adiabatic-like collapse the top down mechanism of structures 
formation is suppressed as soon as viscous effects are taken into account. Such a feature is not present 
in the isothermal-like collapse because the sub-structures formation is yet present and outlines the 
same behavior as in the non-viscous case. We emphasize that in the adiabatic-like collapse the 
bulk viscosity is also responsible for the appearance of a threshold scale beyond which perturbations 
begin to increase; this issue, absent in the non- viscous case, is equivalent to deal with a Jeans length. 
A discussion of the physical character that the choice v = b/& takes place in the present case is 
provided. 

PACS numbers; 95.30.Wi, 51.20.-l-d 



I. INTRODUCTION 



One of the most attractive challenge of Relativistic Cosmology is to provide a self- 
consistent picture for processes of gravitational instability which connect the high isotropy 
of the cosmic microwaves background radiation , 3 with the striking dumpiness of the 
actual Universe An interesting framework is provided by the top-down scheme of struc- 
ture fragmentation, mainly associated with the hot dark matter phenomenology [1], @, 
01 ■ Such a scheme is based on the idea that perturbation scales, contained within a col- 
lapsing gas cloud, start to collapse (forming sub-structures) because their mass overcomes 
the decreasing Jeans value [1], [1], [i^Ij [HI of the background system. The resulting ef- 
fect of such a gravitational instability consists of a progressive enhancement of the density 
contrasts associated to the perturbations sub-scales. 

In a work by C. Hunter a specific model for a gas cloud fragmentation was addressed 
and the behavior of sub-scales density perturbations, outcoming in the extreme collapse, was 
analytically described. The hypotheses on which this model is based are the homogeneity 
and the spherical symmetry, respectively, of the collapsing cloud that starts at rest its 
fall. Furthermore, it is assumed that pressure forces are negligible in the collapse dynamics 
and therefore a real notion of Jeans mass is not required in this approach. According to 
this scheme, in fH is developed the Lagrangian and Eulerian formulation of the zero and 
first-order dynamics respectively. The results of this analysis show that density contrasts 
grow approaching the singularity inducing a fragmentation process of the basic flow. It is 
outlined that pressure effects do not influence the perturbations behavior if we consider an 
isothermal-like politropic index 7 (i.e. for 1 ^ 7 < 4/3). On the other hand, such effects 
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increase as 7 runs from 4/3 to the adiabatic value 5/3. Just this limit case, for 7 = 5/3, 
represents an exception being characterized by density contrasts which remain constant 
asymptotically to the singularity. 

In this work we investigate how the above picture is modified by including, in the gas 
cloud dynamics, the presence of bulk viscosity effects. In this respect, we generalize the 
Lagrangian evolution by taking into account the force acting on the collapsing shell as a 
result of the negative pressure connected to the presence of bulk viscosity. We construct 
such an extension requiring that the asymptotic dynamics of the collapsing cloud is not 
qualitatively affected by the presence of viscosity. In particular we analytically integrate the 
dynamics in correspondence to the constitutive equation [l^] , > [3 > [Hi for the viscosity 
coefficient C = Co where the exponent is assumed to be = 5/6 {p being the cloud 
unperturbed density and Co a- constant parameter). Then we face the Eulerian motion of 
the inhomogeneous perturbations living within the cloud. The resulting viscous dynamics is 
treated in the asymptotic limit to the singularity. We show that density contrasts behave, in 
the isothermal-like collapse, as in the non- viscous case. On the other hand, the perturbations 
damping increases monotonically as 7 runs from 4/3 to 5/3; in fact, for such adiabatic-like 
case, we see that density contrasts asymptotically vanish and no fragmentation processes 
take place in the cloud when the viscous corrections arc sufficiently large. In particular, 
we observe the appearance of a threshold value for the scale of the collapsing perturbations 
depending on the values taken by the parameters Co and 7 G (4/3, 5/3]; such a viscous effect 
corresponds to deal with an analogous of the Jeans length, above which perturbations are 
able to collapse. However such a threshold value does not ensure the diverging behavior of 
density contrasts which takes place, in turn, only when a second (greater) critical length is 
overcome. 

Since in the extreme collapse it is expected that viscous processes are relevant, our anal- 
ysis suggests that the top-down scheme of structures formation can be deeply influenced 
when non-equilibrium features of the dynamics arise. Indeed bulk viscosity outcomes as a 
phenomenological issue of a thermodynamical non-equilibrium collapse and its presence has 
to be consider important in the non-linear regime of structures formation. According to our 
analysis, if such viscous effects are sufficiently intense, the final system configuration is not 
a fragmented cloud as a cluster of sub-structures but simply a single object (a black hole, 
in the present case, because pressure forces are taken negligible). Furthermore, we discuss 
why the choice v = 5/6 has a physical meaning in the viscous dynamics: we show that 
for 1/ > 5/6 the background evolution would be asymptotically affected by viscosity which 
would acquire a non-perturbative character; on the other hand, for v < 5/6 no modifications 
occur with respect to the dynamics of the non- viscous density contrasts. 

The scheme of the paper is as follow. In Sec. [ll] a review of the Hunter work 1] is 
presented; after having integrated the basic flow motion equation we analyze the asymptotic 
(toward the singularity) evolution of the density perturbations in the cases of an isothermal 
and adiabatic behavior of the gas respectively. In Sec. Illll we discuss how the presence of bulk 
viscosity influences the zero-order gas cloud dynamics which is supposed to be qualitatively 
the same as in the non-dissipative case. The analysis of the Lagrangian motion is developed 
for a particular value of the bulk viscosity power-law exponent v = 5/6. In Sec. IIVI we 
start from the unperturbed solutions of the background field to develop a perturbation 
theory in a Newtonian approach. The Eulerian motion equations are written for first- 
order perturbed quantities and then combined to build an unique second-order differential 
equation which governs the density perturbations evolution. In Sec. |V] we integrate this 
fundamental equation to obtain the asymptotic behavior of the fluctuations for the final 
part of the collapse. Two different cases are analyzed: the isothermal-like one which stands 
for 1 ^ 7 < 4/3 and the adiabatic-like case for 4/3 < 7 ^ 5/3. At the end of this Section 
are pointed out the physical reasons to have a viscous fluid with = 5/6 and the conditions 
for the validity of the Newtonian approximation. In Sec. I VII brief concluding remarks follow 
and possible upgradings for the description of the bulk viscosity effects on the gravitational 
collapse are discussed. 
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II. DESCRIPTION OF A NON- VISCOUS CLOUD FRAGMENTATION 



In this Section we present an hydrodynamical analysis of a spherically symmetric gas 
cloud collapse. This model was firstly proposed by Christopher Hunter in 1962 [l], in his 
work he supposed that the gas cloud becomes unstable with respect to its own gravitation 
and begin to condense. The collapsing cloud is assumed to be the dynamical background on 
which studying, in a Newtonian regime, the evolution of density perturbations generated 
on this basic flow. Such an analysis is suitable for the investigation of the cosmological 
structures formation in the top-down scheme Q, Q since it deals with the sub-structures 
temporal evolution compared with the basic flow of the gravitational collapse. 

In the Hunter model, the unperturbed flow was supposed to be homogeneous, spherically 
symmetric and initially at rest. Furthermore the gravitational forces are assumed to be 
very much greater than the pressure ones, which are therefore neglected in the zero-order 
analysis. In such an approach the gas results to be unstable, since there are no forces which 
can contrast the collapse, and the condensation starts immediately. 

The basic flow is governed by the Lagrangian motion equation of a spherically symmetric 
gas distribution which collapses under the only gravitational action. Assuming that the 
initial density of the cloud is constant in space, the dynamics reads 

d\ _ GM 

where the origin O is taken at the center of the gas, r is the radial distance, G the gravita- 
tional constant and M the mass of the gas inside a sphere of radius r. In what follows, we 
shall suppose that the gas was at a distance a from O in correspondence to the initial instant 
to; this distance o identifies a fluid particle and will be used as a Lagrangian independent 
variable so r = r{a,t). Provided that particles do not pass trough each other, the mass 
M inside a sphere of radius r is not time dependent and is a function of a only; using the 
integral form M{a) = /p dr' Anpr''^ = /g da' A'KpQ{a')a''^ , where p is the gas cloud density 
and po = p{to) the initial one, we get the relation 

pr^£=Poa^ (2) 

A first integration of ([T]) yields the expression of the radial velocity v = dr/dt which reads 

V = -[2GM(l/r- l/a)]5 , (3) 

where we considered the negative solution since the radial unperturbed velocity must be 
negative in order to obtain a collapse. Let us now introduce the paramctrization 

r = a cos^ (3 , (4) 

where /? = /3(i) is a time dependent function such that (3{to) — and /3(0) = tt/2 since we 
choose the origin of time to have i = when r — and to takes negative values. Since we 
assumed po to be uniform, we are now able to integrate equation ^ to get the following 
relation between (3 and t and the expression of the initial time tg: 



/3+isin2/3= l+y^TTpoG, (5) 



to = -V3^/32poG. (6) 

It is more convenient to use an Eulerian representation of the flow field. To this end, we now 
derivate the relation ^ and integrate equation ^ to obtain respectively the unperturbed 
radial velocity v and the basic flow density p. Furthermore, solving the well known Poisson 
equation for the gravitational potential 0, we obtain the unperturbed solutions describing 
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the background motion; all these quantities take the explicit forms 



V = [v, 0, 0] , (7a) 

t; = -2r/3tan/3 , (7b) 

p = Po cos"^ (3 , (7c) 

(j) = -27rpoG (a^ - rVS) cos"^ /3 , (7d) 



where (') denotes the derivate with respect to time and in ([Taj) the non-radial components 
of velocity must vanish since we are considering a spherical symmetry. 

The first-order perturbations to the basic flow (higher orders analysis was made by Hunter 
in two later articles [l^, [13]) are investigated in the Newtonian limit starting from the 
continuity, Euler and Poisson equations [l8|. In such a picture, we are now interested to 
study the effects of the thermostatic pressure on the perturbations evolution. We shall 
therefore consider terms due to the pressure p in the motion equations, which read 



p + V-(pv)=0, (8a) 

V + (v • V) • V = -\/(f> --\7p, (8b) 

P 

W^<j> = 4ttGp . (8c) 



The gas is furthermore assumed to be barotropic, i.e. the pressure depends only by the 
background density p. In this model, zero-order solutions ([7]) are already verified since the 
pressure gradient, in the homogeneity hypothesis, vanishes and the pressure affects only the 
perturbative dynamics. 

Let us now investigate first-order fluctuations around the zero-order solutions, i.e. we 
take the perturbed quantities: (v-|-(5v), {p + Sp), {(f> + S(f>) and also {p + Sp) where p = p[p). 
Substituting these solutions in the Euler equation (j8bp and taking the rot of the final 
expression, we get, linearizing in the perturbed quantities, an equation for the vorticity 
Sw — W X Sv which stands 

Sw = —V X [5w X v] . (9) 

Using spherical coordinates [r{a,t),6, ip] we are able to build the solutions for the three 
components of the vorticity, getting 

(5w = [ / cos"** /3 -I- /i , TOCOS"**/?, n cos~'^ f3] . (10) 

Here l,m,n are arbitrary functions of the new variables [a,9,ip] (the radial coordinate 
transforms like ^) which must satisfy the relation V-(5w = ^ and h is physically irrelevant 
since it represents a static distribution of the 5w first component in the space. 

We are now able to find a solution for the perturbed velocity using the vorticity expression; 
one can always consider a solution of the form 

Sv ^Vcos-^ P + V^ , (11) 

where V are arbitrary functions of the coordinates and for V we assume a restriction 
gives by the relation V • V = 0. 

Let us now write three equations for the perturbed quantities (5v, 6p and dcj). Substituting 
the last expression for the velocity fluctuation in equations (|S]) as perturbed to the first- 
order and eliminating the variable r trough the relation ^ we get the following equations 



(5p-6(5p/3tan/3-|-pocos~*°/3i:»V = , (12a) 

il; + 6(f>+^cos^pSp^0, (12b) 
Po 

D^S(P-4ttGcos^ (3 6p = . (12c) 



^ We remember that in any coordinates system the relation div rot = stands. 
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Here time differentiatio n is tak en at some fixed comoving radial coordinate, Vg is the sound 
speed given by Vg — y^dp/dp and is the Laplace operator as written in our comoving 
coordinates. 

A single second-order differential equation for dp can be obtained from the equations 
set ([12]) by eliminating 6(f) between (|12bp and (|12cp and then eliminating ip between the 
resulting relation and equation (I12a|) . This final equation is as follow 



(cos^° pSp-6 sin /3 cos^ /? /3 - 47rGpo cos'' P6p = Vs^ cos^ (3 D^6p . (13) 



In order to study the temporal evolution of density perturbations, we assume to expand 
6p in plane waves of the form 

<5p(r,t) = 5^,(t)e-*-, (14) 

where 1/fc (with k = |k|) represents the initial length scale of the considered fluctuation. 
We shall now express the thermostatic pressure as a function of the basic flow density using 
the barotropical law 

p = Kp^, (15) 

where k, 7 are constants and 1^7^ 5/3. By this expression we are able to distinguish a 
set of different cases related to different values of the politropic index 7. The asymptotic 
value 7=1 represents an isothermal behavior of the gas cloud and corresponds to a constant 
sound speed Vg] the case 7 = 5/3 describes, instead, an adiabatic behavior and it will be 
valid when changes are taking place so fast that no heat is transferred between elements of 
the gas. We can suppose that intermediate values of 7 will describe intermediate types of 
behavior between the isothermal and adiabatic ones. 

The temporal evolution of density perturbations is governed by (|13p : this equation can 
not be solved in general but we can determine the asymptotic behavior of solutions for the 
flnal part of the collapse as (— i) 0. In this limit, we are able to develop up to the first- 
order the equation ([5]) which, once integrated, gives the time dependence of the parameter 
(3. For small /3 we are able to approximate sin /? « 1 in order to obtain the relation 



cos/33 = ^evrGpo {-t) . (16) 

In this approach we determine the asymptotic temporal evolution of the basic fiow unper- 
turbed density (I7cp which now reads 



p-i-t)-'. (17) 

Substituting this expression in (|13p . together with (fH|) and we get the following asymp- 
totic equation for the final part of the collapse 



•• 1 fi 



■ vl P (-t)8/3-27 



5q = Q, (18) 



where Vq = K7P0 ^- ^ complete solution of this equation involves Bessel functions and 
reads 

5b - (-t)-i3/6 [ci J„ [g(-i)'/'"1 + C2 r„ [g(-t)^/3-^]] , (19) 
where the parameters n and q are 

71 = 5/6(4/3-7), (20a) 
q = -«ofc(67rpoG)i/6-''/2 / (4/3 - 7) ■ (20b) 

In order to study the asymptotic evolution of this solution, we shall analyze the cases 1 ^ 
7 < 4/3 and 4/3 < 7 ^ 5/3 separately, since Bessel functions have different limits connected 
to the magnitude of their argument [l9j . In the asymptotic limit to the singularity, the 
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isothermal-like case is characterized by a positive time exponent inside Bessel functions so 
g^4/3-7 ^ 1 on the other hand, in the adiabatic-like behavior we obtain qt^^^^'^ '3> 1. 

For 1 ^ 7 < 4/3 Bessel functions J and Y behave like a power-law of the form J„(x) ~ x+", 
Yn{x) ~ a;-", for a: ^ 1. By this approximation, the solution assumes the the following 
asymptotic form 

SglSO ^ (_^)-3 (21) 

which holds for all the isothermal-like 7 values. This result implies that density pertur- 
bations grow to infinity as (—t) —^ 0; this behavior can be deducted by simply analyzing 
equation p3)) . Its right hand side contribution, due to pressure forces, is proportionally to 
cos^ P which decreases drastically in the last part of the collapse (i.e. (3 7r/2), in this 
way pressure forces become negligible toward the condensation in comparison with respect 
to the gravitational ones and can not, in any cases, prevent the collapse of density fluctu- 
ations. Let us now study the asymptotic behavior of the density contrasts S = 6g/p. It is 
immediate to see that for all the values of 7 inside the interval [1,4/3) density contrasts 
asymptotically diverge like 
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i-t)-' , (22) 



implying that perturbations grow more rapidly with respect to the back-ground density 
favoring the fragmentation of the basic structure independently on the value of the politropic 
index. 

In the adiabatic-like case, for 4/3 < 7 ^ 5/3, the argument of Bessel functions becomes 
much gather than unity and they assume an oscillating behavior like Jn{x) ~ x~^/^ cos{x), 
Yn(x) ~ x~i/2sin(x), for a; ^ 1. The solution (fT9|) asymptotically reads 



Sg^""^ ~ (_^)7/2-17/6 cos 



_(4/3-7)(67rGpo)^/2-i/6 



(23) 



and therefore perturbations oscillate with ever increasing frequency and amplitude. In this 
case, density contrasts assume the form 



:ADB 



i-t)i-e , (24) 



and they outline that perturbations, for intermediate stages as 4/3 < 7 < 5/3, collapse 
before that the basic flow completes the condensation (i.e. 7/2 — 5/6 < 0) and the frag- 
mentation of the background fluid is favored. On the other hand, if the gas cloud behaves 
adiabatically (i.e. 7 = 5/3), perturbations remain of the same order as the basic flow density 
pT]) . We can conclude that, in this adiabatic-like case, pressure forces become progressively 
strong during the collapse as 7 increases having a stabilizing effect which prevents that den- 
sity perturbations grow in amplitude with respect to the unperturbed flow. An intermediate 
type of behavior exists for 7 — 4/3, in this case the disturbances grow like (— t)^^'^/^. 



III. MOTION EQUATIONS OF AN UNPERTURBED VISCOUS FLUID 

In this Section we discuss a model in order to build the motion equations of a spherically 
symmetric and uniform gas cloud including the corrections due to the presence of dissipa- 
tive processes; the hypothesis that the fluid is initially at rest already stands here. The 
Lagrangian equation ([1]) describes a spherical shell which collapses under the gravitational 
action. In such an approach the shell results comoving with the collapsing background, this 
implies that there are no displacements between parts of fluid with respect to ones other 
since we assume an homogeneous and isotropic flow. 

Dissipative processes are therefore related to the only fluid compression and can be phe- 
nomenologically described |18l] by the presence of bulk viscosity (or second viscosity) effects 
summarized by the corresponding coefficient C. Furthermore, in this model we are able to 
neglect the shear viscosity (first viscosity) since it is connected with processes of relative 
motion among different parts of the fluid. 
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The effect of bulk viscosity on the fluid motion can be expressed by the generation of a 
negative pressure additional to the thermostatic one. In the relativistic limit, dissipative 
precesses are derived [l^, [1^, [HI, [l^l as a correction of the energy- momentum tensor 
pressure characterizing the matter filling the space-time. In presence of viscosity, the ther- 
mostatic pressure p is replaced by a new quantity of the form 

p = p~Cu^^, (25) 

where is the shell comoving 4- velocity, i.e. in the present case = (1, 0). The coefficient 
( is not constant and we have to express its dependence on the state parameters of the fluid. 
In the homogeneous model this quantity depends only on time and therefore we may consider 
it as a function of the fluid density p. According to literature developments [T^l, [3) 
[m we assume that ( depends on p via a power-law of the form 

c = Cop^ (26) 

where i^o is a constant and v is an adimensional parameter. 

Let us now evaluate equation (|25p in correspondence to the Newtonian limit appropriate 
to our analysis. As in Sec. |ll] , we consider here that gravitational forces are much greater 
than ones due to the thermostatic pressure in order to neglect p in the zero-order analysis. 
In the non-relativistic approach, the metric we consider is a flat Minkowskian one expressed 
in spherical coordinates, i.e. 

= dt^ - dr^ - r'^dO'^ + sin^ 9 dip^ , (27) 

and the metric determinant is g = — sin^ 9. For the 4-divergence presents in (P5|) we 
immediately obtain, in this case, itf"^ — 2f /r . 

Considering the basic flow density as p = M/{^'nr^), the viscous pressure p now reads 

p = -2Co m/AnY ^ . (28) 

ot 

The pressure force acting on the collapsing shell (and damping its motion) takes the form 
Fp = pAirr'^. By these considerations, the motion equation ([1]) for a viscous fluid becomes 

d^r _ GM_ C dr 

where C = SttCo {iM/Anf. 

The above equation must be integrated to obtain the evolution of the radial velocity v 
and the density p of the unperturbed flow. In order to compare our viscous analysis with 
the Hunter case discussed in the previous Section, let us now require that the viscosity does 
not influence the final form of the velocity and, for instance, it should be yet proportional 
to (see (jSj). Substituting an expression of the form 

V = B/y/? (30) 

into the cc^uation (|29p we see that, in correspondence to the choice v = 5/6, it is again a 
solution as soon as we take the following identification 



B^C- y/C^ + 2GM , i^ = 5/6, (31) 

where B assumes only negative values. Although this dynamics is analytically integrable 
only for the particular value ly = 5/6, the obtained behavior v ^ r~^l'^ remains asymptoti- 
cally (as r ^ 0) valid if the condition < 5/6 is satisfied. 

Using such a solution we are able to build an explicit form of the quantity /3 defined by 
(jH); differentiating this relation with respect to time and taking into account (|30p we obtain 
a differential equation for the variable /3 which admits the solution 

cos = 3A (-i) , (32) 

where A is defined to be A = —Bjlc?!"^. The Eulerian expressions ([7]) of the unperturbed 
quantities v, p and hold here since they are derived simply from relations (|4]) and ([2|): 
the effects of hulk viscosity in this zero order analysis are now only summarized by the new 
time dependence p2p of the parameter (3 which implies a different dynamics of the basic 
flow. 
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IV. PERTURBATIONS THEORY AND THE DYNAMICS OF FLUCTUATIONS 



In the last Section we have discussed the motion of a viscous basic flow which collapses 
under the action of its own gravitation. We shall now suppose that small disturbances 
appear on this field; the perturbations evolution can be described by the Eulerian equations 
dH) modified by the presence of dissipative processes (as in Sec. |ll]we consider here that 
the thermostatic pressure influences the first-order perturbations dynamics). Such kind of 
effects come out from the thermodynamical irreversibility of the collapse process and are 
due to the miscrophysics of non-equilibrium (isj . 

Dissipation effects do not influence the continuity equation (jSap since it is built starting 
from the mass conservation law of the fluid; the Euler equation, on the other hand, requires 
to be modified by the presence of viscosity which causes irreversible transfer of momentum 
from points where the velocity is large to those ones where it is small. A derivation of 



these corrections is developed in [18| starting from general considerations about mechanical 
properties of the fluid. The Euler equation in presence of viscosity corresponds to the 
Navier-Stokes equation when a gravitational potential (j) is included, i.e. 

v-l- (v- V) • V = -V0- i Vp+ - V(V • v) , (33) 

P P 

where C is the bulk viscosity coefficient (|26p (in this formula we have neglected the shear 
viscosity since the homogeneity hypothesis already holds). Beside this equation, the viscous 
fluid motion is described by the continuity equation (|5a|) and the Poisson one ([5c)l . 

Once built motion equations, we start from the unperturbed solutions ([7]) (which still 
hold here as seen in Sec. IIII[) to analyze the behavior, in presence of viscosity, of first-order 
perturbations to the basic flow. Let us now substitute expressions (v-|-5v), {p+5p), {(j)+6(l>) 
and p + Sp into the viscous Euler equation ([55)1 to obtain 

6v + V(v • ^v) + (V X 5v) X V = -VS(f> - -^\/5p + ^ V(V • 6v) . (34) 

P P 

Proceeding as in the non-viscous case (Sec. HII we shall now apply the rot operator to this 
last equation in order to get first the solution for the vorticity Sw and then the expression 
for the velocity perturbations Sv. Indeed, using the vectorial identity rot[V/] — (which 
holds for each scalar function /), all terms in the right hand side of (|34p vanish under 
this operation. In particular the term due to the viscous correction disappears from this 
equation because C is, by assumption, a space-independent function. In this way we reach 
the equation ^ for the vorticity which yields the same solution pT|) as in the non-viscous 
case. 

Following the line of the Hunter work we now build the equations for the perturbed 
quantities (5v, Sp and S<j). Substituting the expression pT|) into the first-order perturbed 
Eulerian motion ([M)) . we obtain (using ^ and the conformal spherical coordinates [a,9,(p]), 
the equation 

iP + S4>+ —cos^ PSp- —cos^ D'^il: = , (35) 
Po Po 

which corresponds to the viscous generalization of (jl2bp . The other perturbed equations 
maintain their own forms (|12aP and (|12cp also in the viscous case. 

Our analysis proceeds in order to build an unique equation which describes the evolution 
of the density perturbations. By following the procedure developed in the non-dissipative 
approach we get now the equation 



( cos^° PSp-6 sin (3 cos^ PPSp^ - AirGpo cos"* /? 6p - 



— sin/3cos"/3/3) D'^Sp+— co?}'^ (3 5p , 
Po > Po 



(36) 



that is the second-order dynamics which governs the density perturbations in the viscous 
regime. 
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V. THE ASYMPTOTIC BEHAVIOR OF DENSITY PERTURBATIONS 



Following the non- viscous approach, let us now factorize perturbations 5p in plane waves 
by the formula (jl4p and then use the barotropic relation p = Kp^ , in which the limit cases 
7=1 and 7 = 5/3 correspond respectively to an isothermal and an adiabatic behavior 
of the gas. According to these assumptions we are able to write the asymptotic form of 
equation ([55)1 near the end of the collapse as <) O.In this case, the quantity cos/3 is 
given by ([32]) . i.e. 



cos/3= (3A)i/3 i-ty/^, 



(37) 



and asymptotically we can make the approximation sin /? « 1 in order to obtain an equation 
which generalizes (fT8|) in presence of viscosity; furthermore substituting expression ([7c|) of 
the basic flow density p into (pS)) . we obtain 



With this assumptions, equation now reads 



(38) 



"14 AnGpo 



16 A 

y y-^) - (3^)2.-2/3 



(-0 



8/3-2iy 



(5^3 - 



+ 



f2 fc2 



(3^1)2^-2/3 



(-0 



8/3-27 



2A 



(3^)2^-2/3 



(39) 



Sg = 



where Vq = ajp^ ^ and the parameter A is given by 

A = Copr'k' 



(40) 



The background motion equations were derived in Sec. IIIII for a particular value of the 
bulk viscosity parameter v — 5/6. In this case, equation (|39p assumes the form 



16 

y 



i-t)dg + 



14 

y 



4ttGpo f2A:2(_<)8/3-27 



9^2 



27-2/3 



2A 
3A 



6g = . (41) 



In analogy with the non-dissipative case, a complete solution of (|4ip involves Bessel 
functions of first and second species J and Y respectively and it explicitly reads 

Sg= CiGiit) + C2G2{t) , (42) 

where Ci , C2 are integration constants and the functions Gi and G2 are defined to be 

G^it) = {~t)-'i + ^ j„[5(_i)4/3-7] ^ G2(i) = {~t)-'i + ^ r„[g(-i)4/3-7] , (43) 

having set the Bessel parameters n and q as 

n = [A^ - 2Ayl + A^ + lenGpo] ^ / (6A(4/3 - 7)) , (44a) 
q = -fcwo(3A)i/^-'^ / (4/3 - 7) . (44b) 

We now proceed, in order to study the asymptotic evolution of the solution (|^^ . analyzing 
the cases 1 ^ 7 < 4/3 and 4/3 < 7 ^ 5/3 separately using the asymptotic expansion for 
Bessel functions introduced in Sec. [Til 

In the first case of an isothermal-like behavior of the gas, the time exponent inside Bessel 
functions results to be positive thus their argument is much less than unity. In this way, an 
asymptotic form of functions G can be found as follow 



Gf° = ci i-ty 



- + S^ + (|-7)r 



Gf« = C2 i-t) 



_n-^+67r-(f-7)r 



(45) 



where Ci and C2 are constants quantities. The condition which implies the density pertur- 
bations collapse is that at least one of G functions diverges as (— t) 0. An analysis of 
time exponents yields that Gi diverges if A < 7 A — 2-kGpq/?>A but, on the other hand, G2 
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is always divergent for all A. These results imply that, in the isothermal case, perturbations 
always condense. Let us now compare these collapses with the basic flow one; the back- 
ground density evolves like cos~^ /? (see (|7c|) ) that is, using ([32|) . in the same way of the 
non-viscous case, i.e. 

p-i-t)-'- (46) 

In the non-dissipative case we have seen that perturbations grow more rapidly with respect 
to the background density involving the fragmentation of the basic flow independently on 
the value of 7; in presence of viscosity density contrasts assume the asymptotic form 

S^SO _ i^__i_yl + ^-J^^[A^-2\A+X'- + W7^Gpo] ^ (47) 

Here the exponent is always negative and it does not depend on 7, this implies that S^^'^ 
diverges as the singularity is approached and real sub-structures are formed involving the 
basic flow fragmentation. This issue means that the viscous forces do not have enough 
strength to contrast an isothermal perturbations collapse implying a formation of an unique 
structure. 

For 4/3 < 7 ^ 5/3 the argument of Bessel functions becomes much gather than unity 
and J, Y assume an oscillating behavior. In this adiabatic-like approach functions G read 

G^^""" - £1,2 {-t)i-'f+^ , (48) 

where C1.2 are constants. Following the isothermal approach, we shall now analyze the time 
power-law exponent in order to determine the collapse conditions. G functions diverge, 
involving perturbations condensation, if the parameter A is less than a threshold value: this 
condition reads A < 17 A — 3A7 for a given value of the barotropic index 7. Expressing A 
in function of the wave number using expression (j40p with = 5/6, we outline, for a fixed 
viscous parameter (^q, a constraint on k which is similar to the condition appearing in the 
Jeans model [^,0]. The threshold value for the wave number is 

kc - ^Jil7A-3jA)pl/'/Co (49) 

and therefore the condition for the density perturbations collapse, i.e. Sp"^^^ — > 00, reads 
k < kc, recalling that, in the Jeans model for a static background, the condition for 
the collapse k < kj = [iirGpo / v-^]^ holds. It is to be remarked that, in absence of 
viscosity (Co = 0), the expression diverges implying that all perturbations scales can be 
conducted to the collapse as in the non-dissipative adiabatic-like case. On the other hand, if 
we consider perturbations of fixed wave number, they asymptotically decrease as (— t) 
for A > 17 A — 3^7. Thus for each k there is a value of the bulk viscosity coefficient over 
which the dissipative forces contrast the formation of sub-structures. 

If fc < kc, perturbations oscillate with ever increasing frequency and amplitude. For a 
non-zero viscosity coefficient, density contrasts evolve like 

5ADB ^ + . (50) 

A study of the time exponent yields the introductions of a new threshold value. If A < 
5A — 3A7, i.e. the viscosity is enough small, sub-structures form; on the other hand, 
when the parameter Co, or the wave number k, provides a A-term overcoming this value, 
the perturbations collapse is so much contrasted that no fragmentation process occurs. In 
other words, if A > 5 A — SA'f we get S"^^^ 0, i.e. for a given 7 there is a viscous 
coefficient Co enough large ables to prevents the sub-structures formation. It is remarkable 
that in the adiabatic case, 7 = 5/3, dissipative processes, of any magnitude order, contrast 
the fragmentation because, while the Jeans-like length survives, the threshold value for sub- 
structures formation approaches infinity. We can conclude that, in the adiabatic-like case 
for 4/3 < 7 ^ 5/3, the fragmentation in the top-down scheme is deeply unfavored by the 
presence of bulk viscosity which strongly contrasts the density perturbations collapse. 

To complete this Section we point out two relevant questions about the generality and 
the validity of our zero and first-order analysis. 
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(i) Physical meaning of v = 5/6 — We now clarify why the choice = 5/6 is appropriate 
to a consistent treatment of the asymptotic viscous collapse. We start by observing that bulk 
viscous effects can be treated in a predictive way only if they behave as small corrections 
to the thermodynamical system. To this respect we have to require that the asymptotic 
collapse is yet appropriately described by the non viscous background flow. As soon as we 
recognize that equation (I29p can be rewritten as follows 

dv C GM 

+ ^ (51) 

it is easy to infer that in the asymptotic limit as r — > 0, the non viscous behavior v ^ r^^l"^ 
is preserved only if 5/6; in fact, in correspondence to this restriction, the viscous correc- 
tion, behaving like Oir^'^^^^l'^)^ is negligible with respect to the leading order 0{r~'^) when 
the singularity is approached; therefore the request that the viscosity is a small correction 
implies the choice v ^ 5/6. 

On the other hand, if we take = 5/6 — A, with A > 0, then the perturbations dynamics 
in the viscous case (P^)) rewrites as 



1^ _ ^ ( fl2A 

3 (3A)i-2A^ ) 



14 _ AnGpo v^P ,n8/3~27 _ 2A 

3 9A2 (3v4)27-2/3 >^ > (3^)1-2A^ ) 



(52) 



Sg^O . 



As far as we deal with the adiabatic-like case, it is immediate to verify that, as {~t) — > 0, 
the viscous terms in 6g and in 6g respectively are negligible and the dynamics matches 
asymptotically the non- viscous Hunter result (|24p (apart from non-relevant features). In 
the isothermal-like case the viscous term in 6g is again negligible, but for A < 4/3 — 7 the 
one appearing in Sg could now dominate the pressure one; however as (— t) both these 
terms provide higher order corrections with respect to the constants in dg and equation 
([5^ reduces to an equation whose solution overlaps the Hunter behavior (we remark 
that for 1/ < 5/6 the viscous parameters asymptotically disappears from the background 
dynamics too). 

Matching together the above considerations for the zero and first-order re-spectively, we 
see that v = 5/6 is the only physical value which does not affect the background dynamics 
but makes important the viscous corrections in the asymptotic behavior of the density 
contrasts. 



(ii) Validity of the Newtonian approximation — Since our analysis addresses Newtonian 
dynamics while the cloud approaches the extreme collapse, it is relevant to precise the con- 
ditions which ensure the validity of such a scheme. The request that the shell corresponding 
to the radial coordinate r lives in the Newtonian paradigm leads to impose that it remains 
greater than its own Schwarzschild Radius, i.e. 

r(t) > 2GM{a) , (53) 

where M{a) = po (|7ra^) and by (|4]) together with the solution p2|) we reache the inequality 

i-t) » -1 (^ttGpo)'^' [s^CoPo - ^{SirCoP^.r + lTrGpoa^-^'^ ] ~' a^/^-a. (54) 

Once fixed the fundamental parameters a, po and the above constraint on the time 
variable states up to which limit a shell remains appropriately described by the Newtonian 
approach. 

About the dynamics of a physical perturbations scale I — (27r/fc) cos/3^ (here 
plays the same role of a cosmic scale factor), its Newtonian evolution is ensured by the 
linear behavior, as soon as, condition ((54)) for the background holds. More precisely a 
perturbations scale is Newtonian if its size is much smaller than the typical space-time 
curvature length, but for a weak gravitational field this requirement must have no physical 
relevance. To explicit such a condition, we require that the physical perturbations scale 
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is much greater than its own Schwarzschild Radius, which leads to the inequahty k ^ 

1/2 

x(— i)^^/'^, where x — [f (27'')^G'/9o(3A)~^/'^] ; combining this result with the inequality 
()54p we arrive to the following constraint 



k » 



277 



(3A)8/3 



(55) 



The condition (l55l) tells us 



being = -1 a3''-3/2 ^g^^^^. _ ^(8^(op-)2 + s^Gp^ ^3-6. 
which modes are Newtonian within the shell whose initial radius takes the value 



VI. CONCLUDING REMARKS 



Our analysis outlined how the presence of bulk viscosity induces a deep modification of 
the extreme gravitational collapse relative to an uniform, spherically symmetric and dust- 
like gas cloud. While the isothermal-like collapse (1 ^ 7 < 4/3) is characterized by sub- 
structures formation even when viscous effects are taken into account, the adiabatic-like one 
(4/3 < 7 ^ 5/3) undergoes an opposite asymptotic regime as soon as the viscosity become 
sufficiently intense. Though bulk viscosity does not affect (by hypothesis) the extreme 
collapse of the background flow, nevertheless its presence changes drastically the dynamics 
of perturbations which are damped at the point to generate vanishing density contrasts. 
Thus, in the adiabatic case, the fate of a collapsing cloud is sensitive to the viscous effects 
by itself induced. In particular bulk viscosity is able to restore a kind of Jeans length for the 
cloud perturbations; scales above this threshold begin to collapse but, if below the second 
threshold, no sub-structures formation takes place. 

The interest in the model here presented comes out because it is expectable that extreme 
regimes of collapse are associated with bulk viscosity which is a macroscopic effect of the non- 
equilibrium thermodynamics characterizing microphysical precesses. However this study 
has to be upgraded by including relativistic effects which becomes asymptotically relevant 
for the small scales. The idea is to restate the problem in terms of a system dynamics 
based on the Einstein equations adapted to the present context. In this respect it would 
be worth to consider radial inhomogeneities of the gas cloud and to include thermostatic 
pressure contributions already on the background flow evolution. Such a generalized scheme 
would allow to evaluate how important bulk viscosity is in favoring black holes or compact 
astrophysical objects with respect to the fragmentation process in sub-structures. We will 
address such line of thinking, relevant for cosmological and astrophysical topics, as subject 
for further investigations in the bulk viscosity phenomenology. 
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